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We present an exact diagonalization study of the single particle spectral function in the ID and 
2D t—J model. By studying the scaling properties with J and t we find a simple building pattern in 
ID and show that every spectral feature can be uniquely assigned by a spinon and holon momentum. 
We find two types of low energy excitations: a band with energy scale J and high spectral weight 
disperses upwards in the interior part of the Brillouin zone and reaches Ef at kp, and a band with 
energy scale t and low spectral weight disperses downwards in the outer part of the zone, touching 
Ef at 3feF . An analogous analysis of the 2D case at half filling shows that the i-band exists also in 
this case, but is diffuse and never reaches the Fermi energy. For the doped case in 2D the picture is 
more reminiscent of ID, in particular the 'main-band' with a dispersion oc J and the 'shadow band' 
with energy scale t can be identified also in this case. This leads us to propose that the shadow 
bands discovered by Aebi ei al. in Bi2212 are the 2D analogue of the Sfc^ singularity in ID systems 
and unrelated to antiferromagnetic spin fiuctuations. 
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I. INTRODUCTION 

Recent developments in photoelectron spectroscopy 
have challenged the apparent simple truth that the Fermi 
surface of cuprate superconductors is simply the one cor- 
responding to LDA band structures with the only ef- 
fect of the closeness to the Mott-Hubbard insulator be- 
ing a moderate correlation narrowing of the band width. 
The discovery of the 'shadow bands' ||l|,^, the temper- 
ature dependent pseudogap in the underdoped state ||] 
and the substantial doping dependence of the quasipar- 
ticle band structure M leave little doubt that a simple 
single-particle description is quite fundamentally inade- 
quate for these materials. Moreover, photoemission ex- 
periments on one-dimensional (ID) copper oxides |q| have 
shown very clear signatures of spin charge separation. 
The equally clear nonobservation of these signatures in 
the cuprate superconductors at any doping level advises 
against another apparent simple truth, namely that the 
Fermi surface seen in the cuprates is simply that of the 
'spinons' in a 2D version of the Tomonaga-Luttinger liq- 
uid (TLL) realized in ID. Motivated by these develop- 
ments, we have performed a detailed exact diagonaliza- 
tion study of the electron removal spectrum in the ID 
and 2D t—J model. This model reads 
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There by the 'constrained' Fermion operators are written 
as Ci^cr = Ci,o-'^i,ff and Si denotes the spin operator on site 
i. The summation (i, j) extends over all pairs of nearest 
neighbors in a ID or 2D square lattice. 
The electron removal spectrum is defined as 
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where Eq and |^o) denote the ground state energy and 
wave function. For small finite clusters, this function can 
be evaluated numerically by means of the Lanczos algo- 
ritm 1^. 

In ID the t—J model is solvable by Bethe ansatz in the 
case J—2t M, but even for this limit the complexity of 
the Bethe ansatz equations precludes an evaluation of 
dynamical correlation functions. For the closely related 
Hubbard model in the limit J/t^O the Bethe-ansatz 
equations simplify ||8|, and an actual calculation of the 
spectral function becomes possible [p||lO|. In all other 
cases Lanczos diagonalization is the only way to obtain 
accurate results for A(fc, a;) [ [ll| . 

In order to analyze our numerical results, we first want 
to develop an intuitive picture of the scaling properties 
of the elementary excitations in ID, which will turn out 
to be useful also in 2D. It has been shown by Ogata and 
Shiba gl that for J/t— >0 the wave functions can be con- 
structed as products of a spinless Fermion wave function, 
which depends only on the positions of the holes, and a 
spin wave function, which depends only on the sequence 
of spins. A naive explanation for this remarkable prop- 
erty is the 'decay' of a hole created in a Neel ordered spin 
background into an uncharged spin-like domain wall, and 
a charged spinless domain wall. Then, since it is the ki- 
netic energy ^t which propagates the charge-like domain 
walls, whereas the exchange energy ~ J moves the spin- 
like domain walls, one may expect that the two types of 
domain walls have different energy scales. Namely the 
excitations of the charge part of the wave function (i.e., 
the 'holons') have t as their energy scale, whereas those of 



the spin part (i.e., the 'spinons') have J as their energy 
scale. Scanning the low energy excitation spectrum of 
ID t~J rings then shows that indeed most of the excited 
states have excitation energies of the form a-t + b- J [|2| , 
which indicates the presence of two different elementary 
excitations with different energy scales. 
Surprisingly enough the low energy spectrum of the 2D 
model shows the same scaling behavior of the excitation 
energies as in ID ||l2], which seems to indicate the exis- 
tence of two types of spin and charge excitations if very 
different nature also in this case. Other cluster results 
indicate, however, that these two types of excitations do 
not exist as 'free particles': the dynamical density cor- 
relation function, which corresponds to the 'particle-hole 
excitations' of holons and shows sharp low energy peaks 
in ID 0] is essentially incoherent in 2D and has prac- 
tically no sharp low energy excitations ||l^. The opti- 
cal conductivity in 2D shows an incoherent high energy 
part with energy scale J ||l4[| - which is completely unex- 
pected for the correlation function of the current operator 
which acts only on the charge degrees of freedom. There 
is moreover rather clear numerical evidence [lq-0| that 
the hole-like low energy excitations can be described to 
very good approximation as spin 1/2 'spin bags' [Q - 
i.e., holes dressed heavily by a local cloud of spin excita- 
tions. 

To obtain further information about similarities and dif- 
ferences between ID and 2D, also in comparison to the 
spectroscopic results, we have performed a systematic 
comparison of the electron removal spectra in both cases. 
As will become apparent, there are some similarities, but 
also clear differences. We suggest that the main differ- 
ence between ID and 2D is a strong attractive interaction 
between 'spinon' and 'holon' in 2D, which leads to a band 
of bound states being pulled out of the continuum of free 
spinon and holon states. This band of bound states - 
which are nothing but simple spin 1/2 Fermions corre- 
sponding to the doped holes - then sets the stage for the 
low energy physics of the system, i.e., true spin-charge 
separation as in ID never occurs. 



II. ONE DIMENSION, HALF FILLING 

We begin with a discussion of the ID model at half- 
filling. Figure shows the electron removal spectra for 
the 12-site ring. Let us first consider the left panel, where 
energies are measured in units of J. Then, one can dis- 
tinguish different types of states according to their scal- 
ing behavior with t: there is one 'band' of peaks (con- 
nected by the thin full line) whose energies relative to 
the single-hole ground state at fc=7r/2 remains practi- 
cally unchanged under a variation of t, i.e., these states 
have J as their energy scale. As a remarkable fact, this 
'band' abruptly disappears half-way in the Brillouin zone, 
i.e., there are no peaks whose energy scales with J be- 



yond fc=7r/2. This looks like a half- filled free-electron 
band with a Fermi level crossing at 7r/2, which however 
is quite remarkable because inverse photoemission is not 
possible at half- filling. Next, in addition to this ' J-band', 
there are several groups of peaks whose excitation en- 
ergy shows a very systematic progression with t. Indeed, 
when plotting the same spectra but measuring energies 
in units of t (right panel of Fig. ||) these peaks coalesce, 
i.e., to excellent approximation the energy scale of these 
states is t. This coexistence of states with different en- 
ergy scales can be nicely seen in the 'double peak' for 
t/J—2 and momentum 27r/6: the peak with lower bind- 
ing energy falls into the J-band, the one with the higher 
binding energy belongs to the t-band. The dispersion of 
the i-band resembles a slightly asymmetric parabola with 
minimum near 7r/2 for the low excitation energies that we 
are considering. The states that fall onto this parabola 
correspond to the creation of a spinon with momentum 
kp——T:/2, and a holon of momentum k + 7r/2. Since the 
spinon momentum is fixed, this group of states then sim- 
ply traces out the holon dispersion. On the other hand, 
the ' J-branch' corresponds to the holon momentum be- 
ing fixed at the minimum of the holon dispersion, and 
thus traces out the spinon dispersion. 
This building principle for the spectra can be pushed fur- 
ther. Namely, one might expect that not only kp but any 
spinon momentum may serve as the starting point for a 
complete branch of peaks which trace out the full holon 
dispersion. That this is indeed the case is shown in Fig. 0. 
There, the entire width of the spectra is shown and we 
have chosen the zero of energy at the excitation energy 
of either the topmost ' J-peak' at kp = tt/2 (left panel) 
or the topmost 'J-peak' at tt/3 (right panel). Due to this 
choice of the zero of energy, the energy oc J of the spinon 
with the respective momentum drops out. Then, when 
measuring energies in units of t different holon bands 'be- 
come sharp', i.e., their energy relative to the respective 
spinon energy scales accurately with t. Moreover, these 
different groups of peaks to good approximation all trace 
out the same simple backfolded nearest neighbor hop- 
ping dispersion, i.e., the dispersion of the holon is simply 
2tcos{kx). As discussed above, the first holon band is 
shifted by the spinous' Fermi momentum, kp = 7r/2, so 
that its dispersion near the band minimum at k — t: 
could be seen in Fig. [|. We have also verified that by 
alligning the spinon peaks at 7r/6 yet another complete 
holon band can be identified. 

We can thus infer the following building principle for 
the spectral function: the basis for the whole con- 
struction is the 'half-filled' spinon band, with dispersion 
— 0.65Jcos(fc); this is indicated by the thick dashed line 
in Fig. Ha. Then, each (fc,(jj)-point of this band provides 
the 'basis' for a complete holon band 2icos(fc), which is 
'hooked on' to the spinon band at its band maximum; 
these holon bands are indicated by the thin full lines in 
Fig. 0(a). Comparison with the numerical results (in this 



case for the 20-site ring) in Fig. g(b) shows that indeed 
to excellent approximation the poles of the single par- 
ticle spectral function fall onto these bands. There are 
some deviations at high binding energies, which however 
are most probably a deficiency of the Lanczos spectra, 
which are highly accurate only at low excitation energy. 
Moreover, the holon bands in Fig. p(b) have been slightly 
shifted, i.e., they are 'hooked on' to the spinon band not 
precisely at their maximum - we have verified that this 
shift has oscillating sign for different chain lengths, so 
that it probably is a finite size effect. As an interest- 
ing feature, the pole strength seems to be constant along 
each of these holon bands, i.e., the weight is a function 
only of the spinon momentum (this seems not to be cor- 
rect for fc=0 and k^n; here it should be noted that for 
these momenta the holon band intersects itself, which 
leads to a doubling of the peak weight). In the ther- 
modynamic limit, the density of bands increases, while 
simultaneously their spectral weight decreases, resulting 
in incoherent continua. Comparing with the exact results 
of Sorella and Parola [g| for the case J/t—^0 it is obvi- 
ous that the outermost holon band in our calculation, 
originating from the spinon Fermi momentum, develops 
into a cusp-like singularity of the spectral weight. The 
spinon band itself, whose energy scale is J, turns into a 
second dispersionless cusp in this limit, which skims at 
zero excitation energy between fc=0 and k—Tr/2. Sorella 
and Parola found the excitation energy of the dispersive 
cusp to be —2t ± 2isin(fc), which corresponds to a the 
backfolded and shifted nearest neighbor hopping band, 
-2i + 2tcos(fc-7r/2). 

Summarizing the data for ID we see that the entire elec- 
tron removal spectrum obeys a very simple building prin- 
ciple, which moreover holds for all momenta and frequen- 
cies. Analyzing the scaling of the different features with J 
and t one can identify 'branches' of states which trace out 
the dispersion of the true elementary excitations of the 
TLL, namely the collective spin and charge excitations. 
The dispersions of the spinons and holons are both con- 
sistent with simple nearest neighbor hopping bands, the 
spinons moreover have a half-filled Fermi surface. While 
these results may not be really new or surprising, we note 
that they demonstrate that exploiting the scaling proper- 
ties of excitation energies provides a very useful method 
to identify the different 'subbands'. In the following, we 
will make extensive use of this principle to address the 
far less understood problems of 2D and finite doping. 



III. TWO DIMENSIONS, HALF FILLING 

We proceed to the 2D model, and also consider first 
the case of half filling. The spectra shown below refer to 
the standard 20-site cluster, which is the largest cluster 
for which the calculation of the electron removal spec- 
trum is feasible also in the doped case. The fc-net for 



this cluster, which is shown in Fig. ^, consists of the 
group of momenta which roughly follows the (1,1) direc- 
tion, and a second group along (tt, 0)^(0, tt). We would 
like to stress that results for other clusters are completely 
consistent with those for the 20-site cluster. Then, the 
left panel of Fig. H shows the photoemission spectrum 
for this cluster at half-filling; thereby we again focus on 
energies within a few J from the top of the band and 
measure energies in units of J. When the spectra are 
aligned at the top of the band, the positions of the other 
dominant low energy peaks do not show a strong varia- 
tion with t. Some peaks do show a slight but systematic 
drift with t, which however is much weaker than in ID. 
A peculiar feature is the peak at (0,0), whose relative 
excitation energy decreases rather than increases with t. 
Inspection shows, however, that the (very weak) disper- 
sion along the line (tt, 0)^(0, tt) (i.e., the lowest three 
momenta in Fig. g|) scales with t to good approximation. 
A possible explanation is the fact that a hole in a 2D sys- 
tem has two distinct mechanisms for propagation, firstly 
by 'string truncation', which gives effective hopping in- 
tegrals ~ J, and secondly by hopping along spiral paths 
[19 1, which gives (smaller) effective hopping integrals ^t 
[20|. It can be shown [gO[ that the dispersion relation for 
a single hole to good approximation can be written as 

E{k) = J ■ ci{cos{kx) + cos{ky))'^ — t ■ C2 cos{kx) cos{ky) 

where ci ;3> C2 > are numerical constants. The 
first term, which originates from the string truncation 
mechanism, gives a dispersion which is degenerate along 
(tt, 0)^(0, tt) and this degeneracy is lifted by the second 
term which is the contribution from the spiral paths; this 
naturally explains the scaling of the dispersion along this 
line with t. Comparing with ID we note that with the 
exception of {tt,tt) the 'J-band' is present in the entire 
Brillouin zone, i.e., the spinon Fermi surface seen at half- 
filling in ID does not exist. 

We turn to the right panel of Fig. |^(b), which shows the 
entire width of the spectra, with energies measured in 
units of t. It is first of all quite obvious that the spec- 
tra generally are more 'diffuse' than in ID, with sharp 
features existing only in the immediate neighborhood of 
the top of the band (except for one relatively sharp high 
energy peak at (7r,7r)). Next, among the diffuse features 
at high energy there are some whose energy accurately 
scales with t. Although these 'peaks' are rather broad, so 
that the assignment of a dispersion is not really meaning- 
ful, their centers of gravity can be roughly fitted by the 
expression —2t ± 2isin(|fcj:| -I- \ky\), which is reminiscent 
of the dispersion of the 'holon-cusp' found by Sorella and 
Parola Q in ID. An important difference as compared 
to ID is the fact that this t-band docs not seem to reach 
the top of the photoemission spectrum - rather it stays 
an energy of ^t below the J-band, which forms the first 
ionization states. We believe that 'in ID language' the 



most plausible interpretation of the data is the forma- 
tion of bound states of spinon and holon: assuming a 
strong attraction between these two excitations, which 
may originate, e.g., from the well-known string mech- 
anism for hole motion in an antiferromagnet [ pi| , one 
may expect that a band of bound states is pulled out of 
the continuum of free spinous and holons. This band of 
bound states corresponds to the J-band (which however 
has a small contribution oc i in its dispersion due to the 
spiral path mechanism). Such a bound state of spinon 
and holon should be a spin-bag-like spin 1/2 Fermion, 
i.e., a hole heavily dressed by spin excitations. There is 
strong numerical evidence (l5-17|, that this is indeed the 
character of the low energy states in 2D at low doping. 
One may expect, however, that such a bound state may 
not be stable for all momenta, and we believe that this 
is the reason for the absence of a J-peak at (7r,7r). In 
this picture, the 2D analogue of the holon is not a co- 
herently propagating excitation, because it is bound to 
the much slower spinon by the linearly ascending string 
potential. This picture fits nicely with the diffuse char- 
acter of the dynamical density correlation function in 2D 
p3[ : this function, which in a TLL should measure ba- 
sically the response of the free holons, in 2D has almost 
exclusively diffuse high energy 'peaks', with virtually no 
sharp low energy peaks. Moreover, the unexpected (in 
the framework of spin charge separation) appearance of 
J as energy scale in the optical conductivity is also readily 
understood in terms of the dipole-excitations of a bound 
spinon- holon pair |lj] . 

Summarizing the data for 2D, we see a band of quasi- 
particle peaks, which predominantly has J as its energy 
scale, and some diffuse high energy 'band' with energy 
scale t. Both, the absence of the 'spinon Fermi surface', 
as well as the lack of sharp 'holon bands' are in clear 
contrast to the situation in ID. The formation of bound 
states of spinon and holon, resulting in a split-off band of 
spin-bag-like spin 1/2 Fermions explains this in a natural 
way. 



IV. ONE DIMENSION, DOPED CASE 

We return to ID and consider the doped case. Figure 
g shows the spectral function for the 12-site ring with 2 
holes. Measuring excitation energies in units of J (left 
panel) we can again identify the spinon band. For 2 holes 
in 12 sites the nominal Fermi momentum is kp=b-K /12 
(i.e., halfway between 7r/3 and 7r/2) and the spinon band 
extends up to this momentum. As was the case at half- 
filling, some other peaks show a systematic progression of 
their excitation energy, and switching the unit of energy 
to t (right panel) again makes a nearly complete 'holon 
band' visible to which these peaks belong. The holon 
band again takes the form of a backfolded tight-binding 
band, but this time the top of the parabola around fc=7r/2 



is missing. The holon band now seems to touch the Fermi 
energy at kp and at 3fcir=97r/12 (the latter momentum 
is half way between 47r/6 and 57r/6). This picture of 
the spectral function nicely fits with the recent exact 
calculation in the limit J/t— >0 by Penc et al. llQ|: on 
the photoemission side, this calculation showed a high 
intensity 'band' which is very similar to the backfolded 
tight-binding dispersion of the holon band. In addition 
there was a dispersionless low intensity band at zero ex- 
citation energy, which corresponds to the spinon band in 
the limit J/t^O. For both, the exact result in the limit 
J/f^O+, and our numerical data for finite J, there are 
thus two branches of states which reach excitation en- 
ergy zero: the 'main band' which touches Ep at kp, and 
the 'shadow band' |1^], which reaches Ep at 3kp. The 
'Fermi level crossings' of these two bands may be thought 
of producing the well known (marginal) singularities in 
the electron momentum distribution n(fc) at kp and 3kp, 
found by Ogata and Shiba |8|. 

The numerical spectra demonstrate a peculiar feature of 
the TLL, namely a kind of Pauli exclusion principle which 
holds for both holons and spinous: the dispersions of both 
types of excitations become incomplete upon doping, i.e., 
the spinon Fermi surface shrinks as if the spinous were 
spin 1/2 particles, while simultaneously the top of the 
holon band is 'sawed off' as if the holons were spinless 
Fermions. It should be noted that this is quite naturally 
to be expected in that the rapidities for the different 'par- 
ticles' in the Bethc ansatz solution both obey a Pauli-like 
exclusion principle M. This has negative implications 
for, e.g., slave boson mean-field calculations, which nec- 
essarily have to treat one type of excitation as a Boson. 
While spin-charge separation is often quoted as justifica- 
tion for the mean-field decoupling, it is obvious that this 
approximation must fail to reproduce the excitation spec- 
trum even qualitatively in ID, the only situation where 
spin-charge separation is really established. 
For a more quantitative discussion of the Fermi points, 
we note that the Fermi momentum for hole concentra- 
tion S is kp — ^{1 — S). For this momentum the first 
branch of low energy excitations reaches Ep. For small 
6 the second branch of low energy excitations comes up 
to Ep at ~3kp -I- 27r = -1(1-1- 3(5). The two marginal 
singularities thus enclose a hole pocket of length 2tt6 
as one would expect for holes corresponding to spinless 
Fermions. It is easy to see, that this hole pocket is noth- 
ing but the manifestation of the holon 'Fermi surface' 
around k=Tr: the lowest charge excitations, which may 
be thought of as corresponding to a particle hole exci- 
tation between the two edges of the holon pocket have 
wave vector Akp — 2tt — 2ttS, i.e., the holon pocket has 
a diameter of 27r/(5, precisely the distance between the 
two marginal singularities. The spectral function for the 
doped case thus follows the same building principle as for 
the case of half filling, with the sole difference being that 
occupied spinon or holon momenta are no longer avail- 



able for the construction of final states. The singularities 
in n{k) may be thought of as enclosing a hole pocket 
corresponding to spinless Fermions, and thus reflect the 
Fermi surface of the holons. The two holon pockets are 
placed such that their inner edges at ±kp enclose the vol- 
ume corresponding to the Fermi sea of spinons of density 
(l/2)(l-<5). 



V. TWO DIMENSIONS, DOPED CASE 

We proceed to the doped case in 2D. Let us note from 
the very beginning that for very simple technical reasons 
the situation is much more unfavorable in this case. To 
begin with, due to the higher symmetry of 2D clusters 
the available k meshes are much coarser: for example, 
amongst the 18 allowed momenta in the 18 site cluster 
only 6 fc-points are actually non symmetry-equivalent, so 
that the amount of nonredundant information is much 
smaller than in ID. Next, unlike ID where a unique rela- 
tionship exists between hole density and Fermi momen- 
tum, most electron numbers in small 2D clusters cor- 
respond to open-shell configurations with highly degen- 
erate ground states for noninteracting particles. In an 
open-shell situation multiplet effects are guaranteed to 
occur, so that it is in general unpredictable which mo- 
menta are occupied and which ones are not (this holds 
for a Fermi liquid, but is most probably true also for 
other 'effective particles'). Unexpected problems may 
arise from this. Bearing this in mind, one therefore may 
not expect to see a similarly detailed and clear picture as 
in ID. 

Then, Figure H shows the photoemission spectra for the 
18-site cluster, with two holes. We first consider the left 
hand panel, where energies are measured in units of J. 
Comparing with Fig. H, some similarities are quite ob- 
vious: the excitation energies of the topmost peaks at 
(0,0), (27r/3,0) are independent of i (although the spec- 
tra for t/J=2 show a slight deviation) so that we can 
identify a 'band' of states with energy scale J. The sit- 
uation actually is not entirely clear, in that the peak at 
(tt/S, 7r/3) is so close in energy to the one at (27r/3,0) 
that it is not possible to decide if their energy differ- 
ence scales with J or t. Next, the topmost peaks at 
(27r/3, 27r/3) and {tt/3, tt) show a systematic progression 
with t, which is very reminiscent of, e.g.. Fig. n^. Plot- 
ting the same spectra with energy scale t indeed to good 
approximation aligns these peaks (although the peak at 
(27r/3, 27r/3) still has a slight drift), i.e., their excitation 
energy relative to the topmost peak at (27r/3,0) scales 
with t. Moreover one can identify a number of diffuse 
'features' at energies between —0.5t and —t, which also 
are roughly aligned; these are indicated by the dashed 
line. In analogy with ID, we can thus distinguish dif- 
ferent branches of states, with different energy scales in 
their excitation energies. While the coarseness of the k- 



meshes introduces some uncertainty, the data are consis- 
tent with a 'J-band' dispersing upwards in the interior 
of the antiferromagnetic Brillouin zone, and a 'i-band' 
dispersing downwards in the outer part, i.e., the same 
situation as seen in ID. A major difference is the fact 
that the 'features' at higher binding energies are all very 
diffuse, at least for J/i>3. More significantly, despite the 
fact that its energy scale seems to be t, the dispersion of 
the 'shadow band' is much weaker than in ID. In other 
words, the effective mass of that band is ~i^^, but with 
a very large prefactor. 

We proceed to the 20-site cluster, also doped with two 
holes (see Fig. g). Choosing J as the unit of energy, we 
see the already familiar situation: the topmost peaks for 
the states at {2tt/5, tt/S) and (tt/S, 37r/5) are aligned (al- 
though t/J=2 again deviates slightly) and several other 
peaks show a systematic progression with t (an unex- 
pected exception is (0, 0) where a well defined peak actu- 
ally is not observed). Changing to energy scale t aligns a 
number of these peaks, which suggests that these peaks 
form a 't-band' which originates from the topmost peak 
at (27r/5,7r/5). This is a second unexpected feature of 
the 20-site cluster, in that for the spectra in ID (and 
for those of the 18-site cluster in 2D) the most intense t- 
band always seemed to originate from the topmost peak 
of the photoemission spectrum. We can only speculate 
that these unusual features are the consequence of, e.g., 
the multiplet effects mentioned above. We also note in 
this context that the spectra at (0, 0) look actually quite 
different for 18 and 20 site cluster, which shows the im- 
pact of finite-size effects. 

Ascribing the special behavior at (0, 0) to finite size- 
effects, we have a quite similar picture as in the 18-site 
cluster, i.e., the topmost peaks for spectra inside the anti- 
ferromagnetic zone have J as their energy scale whereas 
the topmost peaks in the outer part of the zone have 
energy scale t (this also holds for (tt, 0) which is on the 
boundary of the antiferromagnetic zone). As was the case 
in the 18 site cluster the 'shadow band', while having t as 
its energy scale, has a much weaker dispersion than in ID. 
Indeed, fitting the t-bands in both 18 and 20 site cluster 
by an expression of the form 2teff{cos{kx) + cos{ky)) re- 
quires to choose teff « O.lt - it is tempting to speculate 
that this may actually be St, as one would expect e.g. in 
the Gutzwiller picture. Another notable feature is that 
the i-band is restricted to the outer part of the Brillouin 
zone. Only the diffuse high-energy 'band' indicated by 
the dashed line in Figure ^ seems to scale with t. 
For completeness we would like to mention that a similar 
analysis was not possible for the 16-site cluster with 2 
holes. The reason is essentially that for some momenta 
there are no more sharp 'peaks', but rather a multitude 
of densely spaced small peaks. Due to this, we were not 
able to assign any defined 'bands', or groups of peaks 
which showed a systematic scaling of their excitation en- 
ergy. We have also performed this kind of analysis for 



the 16 site cluster with 4 holes and found no more in- 
dication of the energy scale J: at this somewhat higher 
concentration the entire spectra scale with t. 
Summarizing the data for 2D, hole doping seems to lead 
to behavior which is more reminiscent of ID than for half 
filling, in that the J-band dispersing upwards in the in- 
ner part of the Brillouin zone and the i-band dispersing 
downwards in the outer part seem to exist also in this 
case. Much unlike ID, however, the shadow band, while 
in principle having t as its energy scale, still has a very 
weak dispersion, so that the band structure in the doped 
case is practically identical to that in the undoped sys- 
tem [g^. We note however, that the fact that the shadow 
band has t as energy scale has profound implications for 
its explanation: there have been attempts to interpret 
the shadow band in Bi2212 as a 'dynamical replica' of 
the main band, created by scattering of quasiparticles in 
the standard tight-binding band from antiferromagnetic 
spin fluctuations [2^. Experimentally, however, the fact 
that the shadow bands are observed also in the overdoped 
compounds |^, where antiferromagnetism is very weak, 
as well as the fact that they do not seem to become more 
pronounced in the underdoped compounds , where an- 
tiferromagnetism is strong, both suggest otherwise. On 
the theoretical side, we believe that our data very clearly 
rule out this interpretation: both the 'main band' and 
the spin correlation function n^ have J as their relevant 
energy scale, and it would be very hard to understand 
how the energy scale of t for the shadow band should 
emerge from a combination of these two types of ex- 
citations. In fact, the relatively accurate scaling with 
physically very different parameters suggests completely 
different propagation mechanisms for the two types of 
excitations. We therefore believe that the shadow band 
is a separate branch of excitations, probably best com- 
parable to the states which produce the 3kp singularity 
in the ID systems. 



VI. DISCUSSION 

In the previous sections we have investigated the pho- 
toemission spectrum for the one and two dimensional 
t—J model. By studying the parameter dependence of 
the spectra, we could in ID identify 'branches' of states 
which trace out the dispersions of the elementary excita- 
tions of the TLL, the spinous and holons. Both elemen- 
tary excitations have a simple nearest neighbor hopping 
dispersion, but with different band width: that of the 
spinous is ~ J, that of the holons ^t. 
In the doped case there are two groups of states which 
touch the Fermi energy (see Figure O). 'Inside' the non- 
interacting Fermi surface, there is a whole continuum of 
bands dispersing upwards to Ep. The uppermost of these 
bands traces out the spinon dispersion and has J as its 
energy scale, the lowermost band traces out the holon 



dispersion and has t as its energy scale. In the thermo- 
dynamic limit these bands degenerate into 'cusps' and 
merge at Ep ■ In the outer half of the Brillouin zone, there 
are only states which have t as their energy scale. These 
reach the Fermi energy at Skp, giving rise to a second 
Fermi point. While the resolution in k and co available 
in our finite clusters is not sufficient to make statements 
about extreme low energy excitations, the positions of 
the singularities in the electron momentum distribution 
as determined from exact solutions clearly shows that 
both branches of states indeed do touch Ep. The two 
singularities may be thought of enclosing a hole pocket 
of extent 2TrS, which is essentially the image of the holon 
Fermi surface. 

In 2D at half-filling, the situation is quite different: while 
it is still possible to distinguish bands with different scal- 
ing behavior with J and t, the spinon Fermi surface 
present in ID does not exist and the 'holons' seem to cor- 
respond to overdamped resonances rather than sharp ex- 
citations as in ID. We propose that strong attraction be- 
tween spin and charge excitations, most probably due to 
the well-known string mechanism, pulls a band of bound 
states out of the continuum of 'free' holon and spinon 
states. The relevant physics thus is that of spin-bag-like 
spin 1/2 quasiparticles, as suggested by a considerable 
amount of numerical evidence. 

For the doped case in 2D, the situation is less clear and 
actually somewhat ambiguous. The numerical photoe- 
mission spectra show some analogy with ID, in that there 
seems to be a high intensity 'main band' with energy scale 
J dispersing upwards in the inner part of the Brillouin 
zone, and a low intensity 'shadow band' with energy scale 
t dispersing downwards in the outer part of the Brillouin 
zone (see Figure P). In contrast to ID the dispersion of 
the shadow band is much weaker, i.e., while the energy 
scale of the dispersion is t, it has an additional very small 
prefactor (of the order of the hole concentration). More- 
over the t-band seems limited to the outer part of the 
Brillouin zone, i.e. there are no indications for a holon 
band with energy scale t dispersing upwards in the inner 
part of the Brillouin zone. Only in the 18-site cluster a 
diffuse 'band' with energy scale t can be roughly iden- 
tified at higher binding energies. The different energy 
scales of main band and shadow band suggest that these 
are excitations of quite different nature, and in particu- 
lar rule out the explanation that that the shadow band 
is created by scattering from antiferromagnetic spin fluc- 
tuations. 

Turning to experiment, the results for 2D immediately 
suggests a comparison with the data of Aebi et al. 
m. These authors found that in addition to the 'bright' 
part of the band structure, which seems to be consistent 
with the noninteracting one, there is also a low inten- 
sity 'replica', shifted approximately by (tt, tt), which had 
been consistently overlooked in all previous studies. If 
one wants to make a correspondence to the situation for 



the t—J model, one thus should identify this low inten- 
sity part with the i-band dispersing downwards in the 
t—J model. Our data imply that the shadow band should 
have a slightly different dispersion than the main band. 
The limitations of the cluster method probably preclude 
any meaningful quantitative statements, but it might be 
interesting to see if this difference in dispersion can be 
resolved experimentally. 

We conclude by outlining a somewhat speculative sce- 
nario, based on the assumption that the two bands rep- 
resent indeed different excitations, which persist at all 
temperatures and independent of antiferromagnetic cor- 
relations. In this case, the topology in 2D opens an in- 
teresting possibility: whereas in ID the two classes of 
low energy excitations forming the kp and Skp singu- 
larities in n(k) are well separated in k and lu space for 
simple topological reasons, the experimental data of Aebi 
et al. indicate that the main and shadow band inter- 
sect at certain points in the Brillouin zone (see Fig. |l^, 
left panel). Neglecting the small difference in disper- 
sion between main and shadow band, we might therefore 
model the low energy excitation spectrum by the effective 
Hamiltonian 



H, 



QP 



E 
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certain temperature T* then could produce a crossover 
from the hole pockets to the 'large' Fermi surface at T*, 
which picture would nicely reproduce the pseudogap phe- 
nomenology observed H in cuprate superconductors. 
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FIG. 1. Photoemission spectra at half-filling for different 
values of J/t. The thin lines mark the spinon and holon band. 



FIG. 2. Photoemission spectra at half-filling in the 12-site 
ring. For all J/t the energy of the state marked by the black 
dot is taken as the zero of energy. 



FIG. 3. (a) Schematic building principle of the spectral 
function at half-filling: Spinon band (thick dashed lines) , and 
holon bands (thin full lines), (b) Theoretical spectrum over- 
laid with the numerical result for the 20-site ring with J/t— 10. 
The centers of the circles give the excitation energies and their 
diameter the pole strength of the peaks in the electron removal 
spectrum. 



where e(fc) is the dispersion of the main band, Q — 
(tt, tt), and the a and b operators refer to the main and 
shadow band, respectively. Choosing a dispersion of the 
form e{k) = —2t{cos{kx) + cos{ky)) + At' cos{kx) cos{ky) 
this Hamiltonian reproduces the Fermi surface topology 
found by Aebi et al. quite well (see Fig. |l^). However, 
as mentioned above the two branches of excitations inter- 
sect at some points of the Brillouin zone, so that already a 
small mixing between the two bands, which in turn may 
originate from the spinon-holon interaction, has a dra- 
matic effect on the topology of the low energy excitation 
spectrum. Namely adding a term of the form 
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i.e., a hybridization between the two types of bands, even 
relatively small values of A open up a gap around (tt, 0) 
and transform the Fermi surface transformed into a hole 
pocket (see left panel of Fig. ^^. Thereby we fix the 
chemical potential by requiring that the number of a and 
b particles remains unchanged; it is easy to see that the 
area covered by the pockets then equals the hole con- 
centration (5, precisely as it was the case in ID. Thereby 
the Fermi surface has predominant main band charac- 
ter at its inner edge, and shadow band character at the 
outer edge, implying a very different 'visibility' in pho- 
toemission. Finally, it is tempting to speculate that the 
'pseudo gap order parameter' A decreases with increas- 
ing temperature/hole concentration. Its vanishing at a 



FIG. 4. Allowed momenta for the 20-site cluster. 



FIG. 5. Photoemission spectra at half-filling for the 2D 
20-site cluster. For all J/t the zero of energy has been set 
to the energy of the peak marked by the black dot. 

FIG. 6. Photoemission spectra for the 12-site ring with 2 
holes. 



FIG. 7. Photoemission spectra for the 2D 18-site cluster 
with 2 holes. The spectra for momenta marked by an asterisk 
have been multiplied by a factor of 2 for clarity. 



FIG. 8. Photoemission spectra for the 2D 20-site cluster 
with 2 holes. The spectra for momenta marked by an asterisk 
have been multiplied by a factor of 2 for clarity. 



FIG. 9. Schematic representation of the 'band structure' 
in ID and 2D systems as deduced from the diagonalization 
spectra. 



FIG. 10. Left hand panel: Fermi surface for the Hamilto- 
nian (U) with t' /t = 0.3, the full line corresponds to the main 
band, the dashed line to the shadow band. Right hand panel: 
Fermi surface for t' /t^Q.i, A/t=0.2. 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 
[7] 

[8] 
[9] 

[10] 

[11] 
[12] 



P. Aebi et al, Phys. Rev. Lett. 72, 2757 (1995). 

S. La Rosa et al, preprint. 

A. G. Loeser et al, Science 273, 325 (1996). 

D. S. Marshall et al., Phys. Rev. Lett. 76, 4841 (1996). 
C. Kim et al., Phys. Rev. Lett. 77, 4054 (1996). 

E. Dagotto, Rev. Mod. Phys. 66, 763 (1994). 

P.A. Bares and M. Blatter, Phys. Rev. Lett. 64, 2567 

(1990). 

M. Ogata and H. Shiba, Phys. Rev. B 41, 2326 (1990). 

S. Sorella and A. Parola, J. Phys. Cond. Mat. 4, 3589 

(1992). 

K. Penc, K. Hallberg, and H. Shiba, Phys. Rev. Lett. 77, 

1390 (1996). 

J. Favand et al, e-print :ond-mat 9611221^ . 

R. Eder et al, unpublished. 



[13] R. Eder, Y. Ohta, and S. Maekawa, Phys. Rev. Lett. 74, 

5124 (1995). 
[14] R. Eder, P. Wrobel, and Y. Ohta, Phys. Rev. B 54, 

R11034 (1996). 
[15] E. Dagotto and J. R. Schrieffer, Phys. Rev. B 43, 8705 

(1990). 
[16] R. Eder and Y. Ohta, Phys. Rev. B 50, 10043 (1994). 
[17] J. Riera and E. Dagotto, unpublished. 
[18] J. R. Schrieffer, X. G. Wen, and S. C. Zhang, Phys. Rev. 

B 39, 11663 (1989). 
[19] S. A. Trugman, Phys. Rev. B 37, 1597 (1987). 
[20] R. Eder and K. W. Becker, Z. Phys. B. 78, 219 (1990); 

see also M. Vojta and K. W. Becker, Phys. Rev. B 54, 

15483 (1996). 
[21] L. N. Bulaevskii, E. L. Nagaev, and D. I. Khomskii, Sov. 

Phys. JETP 27, 836 (1968). 
[22] R. Eder, Y. Ohta, and T. Shimozato, Phys. Rev. B 50, 

3350 (1994). 
[23] A. Chubukov, Phys. Rev. B 52, R3840 (1995). 




-2.5 -2 -1.5 -1 -0.5 -1.5 

cj/J 



-0.5 

cj/t 








-3 -2 

cj/t 







-3 -2 

cj/t 








cj/t 







^^^^^^ 


§ 


1 


1 


(b) 


^^x^^\^\^ 


\ 




^ 




^^\\» \X 


N 


^ 


\ 


\^ 


5v\v\^ 


0\ 


®\ 


\ 


\7\ 


(sN\>x\ 


o\ 


♦\ 


\ 


3^ 


©kV\\ 




0\ 




^\ 


1 1 1 1 


1 


1 


• 
1 


1 



k 



7T 




■1.5 



-1 -0.5 

cj/J 







-4 -2 

cj/t 








-2.5 -2 -1.5 -1 -0.5 -4 -3 -2 -1 

cj/J cj/t 







■1.5 -1 -0.5 -1 

CJ/J 



-0.5 

cj/t 








(0,0) 



jy.''-jj^^ 



l\ 



'■ ,.^V ^ /■ V \.\ 



(27t/5,7t/5) 




i^'XA-yy ^.'^.^v". A7V Av- 



(tt.O) 



^' V .^A.^>>jg^;^-.^iw4L\ A.'u >- 



fill 



(471/5.277/5)" 



^A 



JL 



^;^-V^-- 



JL 



(371/5,471/5)" 



..^aaJ^A. I L^VV 




)'J\^.L-'^^-L:' 



-1.5 -1 -0.5 -0.5 -0.25 

cj/J cj/t 



•-F 


spinon /T 




\ 


^rr\a\rVx 




shadow' 


O) 

^_ 

CD 
C 
JJ 


* 

' ho on 
.' t 

1 




',ho on 
', f 

% 

\ 
\ 

\ 
\ 
1 

\ 


t 

'ho on' 
(diffuse) 




t 




' 


1D 


\ 
\ 




2D 









^(0,0) 



(7l,7l) 



Figure 9 



(tt.tt) 













